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^ . We report the detection of the primary transit of the extra-solar planet HD 80606 b, thanks to photometric and spectroscopic obser- 

5-H ' vations performed at the Observatoire de Haute-Provence, simultaneously with the CCD camera at the 120-cm telescope and the 

j/j ' SOPHIE spectrograph on the 193-cm telescope. We observed the whole egress of the transit and partially its central part, in both 

, datasets with the same timings. The ingress occurred before sunset so was not observed. The full duration of the transit was between 

9.5 and 17.2 hours. The data allows the planetary radius to be measured (Rp = 0.9 ± 0. 10Rj„p) and other parameters of the system 
to be refined. Radial velocity measurements show the detection of a prograde Rossiter-McLaughlin effect, and provide a hint of a 
spin-orbit misalignment. If confirmed, this misalignment would corroborate the hypothesis that HD 80606 b owes its unusual orbital 
configuration to Kozai migration. HD 80606 b is by far the transiting planet on the longest period detected today. Its unusually small 
radius reinforces the observed relationship between the planet radius and the incident flux received from the star and opens new 
questions for theory. Orbiting a bright star (V = 9), it opens opportunities for numerous follow-up studies. 
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^\ 1. Introduction two telescopes of the Observatoire de Haute-Provence (OHP), 

. France: the CCD camera at the 120-cm telescope, and the 

^ ^^^^^'^l'"'"' P u 7^^^Z^'^ ■ SOPHIE spectrograph at the 1.93-m telescope. This allows us 

^ ELODIE spectrograph by Naef et al. mB- This is a giant detection of the primary transit of HD 80606 b, 

planet (4Mj,p) with a 1 11 -day orbital period on an extremely photometry as well as in spectroscopy through the Rossiter- 
j_] eccentric orbit (e = 0.93). HD 80606 is a member of a binary McLaughlin effect 
Ci system (with HD 80607) with a separation of ~ 1200AU(~ 20" We recall the stellar characteristics of the primary star, which 
on the sky). Wu & Murray (2003) suggest that the present orbit ^^^^ ^^^^ ^^^^^-^ ^^^^ presented here- 

of HD 80606 b results from the combination of the Kozai cy- j^j^ gQ^Q^ ^ ^^^^ ^-^^ ^ ^^^jl^^ measured by 

cles (induced by the distant stellar companion) and tidal friction. Hipparcos of 17+5 mas. A compilation of specti'oscopic data 
Recently, Laughlin et al.^ (2009 ) report detecting a secondary ^^^^ literature gives an effective temperature of 5574 + 50 K, 
transit for HD 80606 b usmg 8 //m Spitzer observations around j 4 45 ^ q 05 ^ ^-^ metallicity of 0.33 dex (Santos 

the periastron passage. This imphes an inchnation of the system y^^^^^- &FischerE005). The stellar mass can be es- 

near / = 90 , and a ~ 15 % probabihty that the planet also shows ^-^^^^^ ^^-^^ isochrones and we get 0.98 + 0. 10 M^. Using the 
primary transits. relationship between luminosity, temperature, gravity, and mass. 

If transits occur, opportunities for detecting them are rare jj^^ ^j^jj^^ luminosity is estimated 0.84 + 0.13 L^. Finally, relat- 
because the orbital period of the planet is almost four months. -^^ ^^^^^ ^-^^^ luminosity and temperature, we derive a radius of 
We managed an observational campaign to attempt detection of 0.98+0.07 R^. These mean values were obtained after several it- 
the transit of HD 80606 b scheduled to happen on Valentine s ^^^^-^^^ ^^^^ ^^^^ ^^^^^ determination. Chromospheric ac- 
Day (14 February 2009). We simultaneously used instruments of ^^-^^ ^^^^ ^„ -^^^^^^^ ^ 7 7 ^ Gyrs (Saffe et al. .20051 

* Based on observations made with the 1 .20-m and 1 .93-m telescopes 
at the Observatoire de Haute-Provence (CNRS), France, by the SOPHIE 2. The phOtOITietriC transit Of HD 80606b 

consortium (program 07A.PNP.CONS). Electronic tables of the data 

are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr Predicted epochs for th e prim ary transit of HD 80606b were 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+ A/, given by Laughlin et al. (120091 ). We carried observations around 
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the expected transit epoch JD=2,454,876.5 with the 120-cm tele- 
scope at OHP, equipped with a 12 arcmin x 12 arcmin CCD 
camera. The Bessel R filter and a neutral density were inserted, 
to insure unsaturated, focused images of the V - 9 target. We 
obtained 326 frames on 13 Feburary 2009 and 238 frames on 
the preceeding night, for comparison. Typical exposure times 
range between 60 sec on the first night and 20-30 sec on the 
second. Aperture photometry was then performed on both data 
sequences. Apertures of 8 and 6 pixels were used for the first and 
second observing nights, respectively. The secondary compan- 
ion HD 80607 is taken as a reference for HD 80606. Both stars 
are separated by 24 pixels, which prevents contamination even 
using simple aperture photometry. The sky background is evalu- 
ated in rings of about 12-15 pixel radii. The resulting lightcurve 
is shown in Fig. [T](upper panel) with all data included. The data 
quality is significantly better during the transit night, because of 
different seeing conditions. The nns is about 0.0023 and 0.0030, 
respectively. 

An egress is clearly detected in the data sequence obtained 
during the night 13-14 February. A shift of almost one half tran- 
sit is observed, in comparison to expected ephemeris. Long-term 
systematics are observed in the lightcurve and removed by a 
polynomial function of the airmass, with the criterion of get- 
ting a flat section of the out-of-transit flux. This correction does 
not strongly affect the transit shape. It is checked on the 12-13 
February sequence that long-term fluctuations are low (not cor- 
rected for in Fig.lTJ. The beginning of the transit sequence unam- 
biguously shows that we do not detect the ingress of the transit. 
The first hour of the sequence indicates a slight decrease but the 
data are quite noisy due to the low object's elevation, and this 
may be introduced by the coiTection for airmass variations. We 
observed in total 7 hours during transit on the second night, and 
3.4 hours after the transit. In addition, we gathered 9.8 hours out 
of transit on the first night. 
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Fig. 1. Photometry (top) and radial velocities (bottom) of 
HD 80606 from 12 to 14 February, 2009, obtained at OHP wifli 
the 120-cm and 193-cm telescopes, respectively. The planetary 
transit is detected in both datasets at the same timing. Top: 
Superimposed are the two extreme {b = in red-dashed, and 
b = 0.91 in green-long-dashed) and the mean (in blue-dotted, 
b = 0.75) models that coiTespond to our data set. Bottom: The 
orbital solution is overplotted (solid line. Table [TJ, together with 
Rossiter-McLaughlin effect models presented in Table |2] in red- 
dashed (^7 = 0, .1 = 0°), blue-dotted (b = 0.75, A. = 63°), and 
green-long-dashed lines (b = 0.91, A = 80°). 



Modelling the primary transit lightcurve of HD 80606 b is 
done in the first place using models of circular orbits, to con- 
strain the inclination and the radius ratio. The Universal Transit 
Modeler (Deeg |20081 l is used, including the limb-darkening co- 
efficients of Claret (2000) for the r' filter, and parameters of the 
orbits given in the next section. Figures [T] and |2] show three 
transit models superimposed to the data, corresponding to im- 
pact parameter b ranging from to 0.91 or inclinations ranging 
from 89.2 to 90.0 deg. The O - C residuals depicted in Fig|2] 
coiTespond to an average transit duration of 13.5 hours, with 
b = 0.75 and nns of 0.0023. The transit depth imposes a ra- 
dius ratio Rp/Rt of 0.094 + 0.009. The planet radius is then esti- 
mated to be 0.9 + 0. 10 Rjup. To match the full transit lightcurve, a 
model including eccentricity would be required. The asymmetry 
of the ingress and egress should be detected and properly fitted, 
for instance. Since we have a partial transit, the approximation 
of the circular modelling is acceptable here, if one takes the rel- 
ative projection of the transit angle and the line of sight into ac- 
count. In a further study, we plan to investigate the modelling of 
the asymetric transit by including the eccentric orbit. We do not 
expect major differences compared to the simple fit performed 
here, before new, more complete photometric data are obtained. 
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Fig. 2. Zoom-in plot on the photometric transit and the three 
models with the impact parameter ranging from to 0.91. The 
residuals are shown below, with an offset of 0.02 on the Y-axis 
for clarity, and they correspond to the mean model with b = 0.75 . 



3. The spectroscopic transit of HD 80606b 

We observed HD 80606 with flie SOPHIE instrument at flie 1 .93- 
m telescope of OHP. SOPHIE is a cross-dispersed, environmen- 
tally stabilized echelle spectrograph dedicated to high-precision 
radial velocity measurements (Bouchy et al. 120061 PeiTuchot et 
al. |2008jl. We used the high-resolution mode (resolution power 
R - 75,000) of the spectrograph and the fast-readout mode of 
the CCD detector. The two 3"-wide circular apertures (optical 
fibers) were used, the first one centred on the target and the sec- 
ond one on the sky to simultaneously measure its background. 
This second aperture, 2' away from the first one, allows us to 
check that there is no significant pollution due to moonlight on 
the spectra of the target. We obtained 48 radial velocity mea- 
surements from 8 to 17 February 2009, including a full sequence 
during the night 13 February (BJD = 54876), when the possible 
transit was expected to occur according to the ephemeris. The 
exposure times range from 600 to 1500 sec, insuring a constant 
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signal-to-noise ratio. This observation was performed in parallel 
to the photometric ones. 

The sequence of the transit night is plotted in Fig. [1] lower 
panel, together with the measurement secured the previous night. 
The Keplerian curve expected from the orbital parameters is 
overplotted. The radial velocities of the 13 February night are 
clearly blue-shifted by ~ 10 ms ' from the Keplerian curve 
in the 1st half of the night, then match the Keplerian curve in 
the 2nd half of the night. This is the feature expected for a 
planet transiting on a prograde orbit, according to the Rossiter- 
McLaughlin (RM) effect. This effect occurs when an object tran- 
sits in front of a rotating star, causing a spectral distortion of the 
stellar lines profile, and thus resulting in a Doppler-shift anomaly 
(Ohta et al. 120051 Gimenez et al. l2006bl Gaudi & Winn|2003l. 

On the RM feature of HD 80606b (Fig. [B, the third and 
fourth contacts occuiTed at BJD ^ 54876.45 and BJD =^ 
54876.55, respectively, whereas the first contact occurred before 
sunset and was not observed. These timings agree with those 
of the photometry (Sect. 2 and Fig. |2]i and the detection of the 
Rossiter-McLaughlin anomaly is unambiguous. 

The Keplerian curve in Fig.[T]coiTesponds to the orbital pa- 
rameters that we refined for HD 80606 b. We used the SOPHIE 
measurements performed out of the transit, as well as 45 Keck 
measurements (Butler et al. 2006) and 74 ELODIE measure- 
ments (55 published by Naef et al. (2001) and 19 additional 
measurements obtained from BJD = 52075 to 52961). We al- 
lowed free radial- velocity shifts between the three datasets. We 
used the constraint of the secondary transit given by Laughlin 
et al. ( [2009] l (T, = 2454424.736 + 0.003 HJD). We also used 
our constraint on the primary transit considering that the end of 
transit is Tigress = 2454 876.55 ± 0.03 BJD. From these con- 
traints, we estimated that the inclination of the system is from 
90° {T, = 2454 876.20 BJD with 17.2 h duration) to 89.2° 
(T, = 2454 876.32 BJD with 9.4 h duration). 

Assuming those constraints, we adjusted the Keplerian orbit. 
The dispersion of the radial velocities around this fit is 8.6 m s ' , 
and the reduced is 1.4. The obtained parameters are reported 
in Table [1] They agree with those of Laughlin et al. ( |2009l l. ex- 
cept for the period, where there is a 3-cr disagreement. The full 
data set and orbital solution are plotted in Fig. [3] We note that no 
anterior data was obtained during the transit by any instrument, 
as shown in the inset of Fig. [3] 

To model the RM effect, we used the analytical approach 
developed by Ohta et al. (2005). The complete model has 12 
parameters: the six standard orbital parameters, the radius ratio 
rp/Rt, the orbital semi-major axis to stellar radius a/R, (con- 
strained by the transit duration), the sky-projected angle between 
the stellar spin axis and the planetary orbital axis A, the sky- 
projected stellar rotational velocity v sin /, the orbital inclination 
i, and the stellar limb-darkening coefficient e. For our purpose, 
we used the orbital parameters and photometric transit param- 
eters as derived previously. We fixed the linear limb-darkening 
coefficient e = 0.78, based on Claret ( |20001 l tables for filter g' 
and for the stellar parameters derived in Sect. 1 . Our free param- 
eters are then A, v sin /, and /. As we observed a partial transit, 
there is no way to put a strong constraint on the inclination /. 
We then decided to adjust A for different values of / in the range 
89.2 - 90 °. The results of our fits (Table |2] and Fig. [H lower 
panel) first show that the stellar rotation is prograde relative to 
the planet orbit. Assuming / = 90° and A - 0°, the projected 
rotation velocity of the star vsin/ determined by our RM fit 
is 2.2 kms"'. This agrees with the value 1.8 kms"' obtained 
by Valenti & Fisher (2005), as well as our spectroscopic de- 
termination (2-3 kms ') from SOPHIE spectra. This latest one 
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Fig. 3. Phase-folded radial velocity measurements of HD 80606 
as a function of the orbital phase, and Keplerian fit to the 
data. ELODIE data in blue. Keck data in red, SOPHIE data in 
green. Orbital parameters coiTesponding to this fit are reported 
in Table[T] The inset shows a zoom around the transit phase. One 
Keck spectrum was obtained Ihr after transit. 



could be slightly overestimated here due to the high metallicity 
of HD 80606. We decided to fix this value and to explore the dif- 
ferent values of inclination angle / to estimate the spin-orbit A 
angle. We see in Table|2]that, if the transit is not central, then the 
RM fit suggests that the spin-orbit angle is not aligned. 



4. Discussion and conclusion 

Despite the low probability of a 111-day period system being 
seen edge-on at both at the primary and the secondary transit 
phases (about 1 % in the case of HD 80606 b), the data acquired 



Table 1. Fitted orbit and planetary parameters for HD 80606b. 



Parameters 


Values and i-cr error bars 


Unit 


Vr (Elodie) 


3.788 ± 0.002 


kms" 


V, (SOPHIE) 


3.911+0.002 


kms" 


P 


111.436 ±0.003 


days 


e 


0.934 ± 0.003 




CO 


300.6 ± 0.4 




K 


472 ±5 


ms"' 


Tq (periastron) 


2454424857 ±0.05 


BJD 


Mp sin i 


40 ± 0.3"^ 


Mj,p 


a 


0.453 ±0.015"^ 


AU 


T, (primary transit) 


2 454 876.27 ± 0.08 


BJD 


(secondary transit) 


2 454 424.736 ±0.003* 


BJD 


tl4 


9.5 - 17.2 


hours 


t23 


8.7-15.7 


hours 




0.98 ±0.10 


Mo 




0.98 ± 0.07 


Ro 


R,/R, 


0.094 ± 0.009 




R, 


0.9 ±0.10 


Rjup 


b 


0.75 (-0.75, +0.16) 


i 


89.6(-0.4, +0.4) 





t: using = 0.98 ± O.lOMj 
p. from Laughlin et al. ( i2009t 
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Table 2. Parameter sets for the Rossiter-McLaughlin effect mod- 
els. 



i 


transit duration 


Ttransit 


spin-orbit A 


x' 


(deg) 


(iiours) 


BJD -2454000 


(deg) 




90.0 


17.2 


876.20 





33.0 


89.6 


15.5 


876.24 


63 


25.3 


89.2 


9.4 


876.32 


86 


33.4 



at the Observatoire de Haute-Provence and presented here show 
this alignment unambiguously. With a partial transit observed, 
we were able to constrain the orbital parameters, including the 
inclination with a precision of ^ 0.4°, and to measure the plane- 
tary radius. The error bars of our measurements should be taken 
with caution, however, and the system's parameters slightly re- 
vised when more complete data is obtained, since systematic 
noise is more difficult to coiTect with incomplete transits. 

The planet has a low radius (0.9 Rjup) considering its mass 
(4 Mjup). Since it is also by far the transiting gas giant receiving 
the lower iiTadiation from its parent star, it is tempting to see its 
small radius as reinforcing the explanation of anomalously large 
hot Jupiters as caused primarily to stellar irradiation, as proposed 
for instance by Guillot & Showman (2002). Figure |4] shows the 
increasingly clear correlation between equilibrium temperature 
and size for transiting gas giants. The relation holds for both 
massive planets (> 2 Mjup) and Jupiter-like planets. Tidal effects 
(Jackson et al., 2007) may play a role in the observed radius of 
HD 80606 b. The high metallicity of the parent star also helps 
provide refractory material for a massive core, although the re- 
quired enrichment would be beyond actual expectations. More 
theoretical development is needed to reproduce the system's pa- 
rameters, taking the whole history of orbital evolution and vari- 
ations in the iiTadiation conditions into account. 

Most of the ~ 60 known transiting planets are orbiting close 
to their hosting stars. Only 5 of them have periods longer than 
5 days, the most distant from its star being HD 17156b, on a 
2 1 .2-day period. HD 80606 b has by far the longest period (111.4 
days). It may be compared to other planetary systems with a 
massive planet in an eccentric orbit: HD 17156, HAT-P-2, and 
XO-3 . HD 80606 b has a smaller radius than those planets, which 
can be related to the migration history or to the changes in 
stellar iiTadiation along the orbit (Laughlin et al, 2009). The 
shape of the Rossiter-McLaughlin anomaly shows that the or- 
bit of HD 80606 b is prograde, and suggests that it could be sig- 
nificantly inclined relative to the stellar equator Since a high 
initial relative inclination is a key requirement for Kozai mi- 
gration to work (Wu & Murray 2003), this observation is not 
surprising. Tighter constraints on the spin-orbit misalignment in 
HD 80606 may support the Wu & Murray formation scenario 
and may provide compelling evidence that the orbital evolution 
of HD 80606 b was once dominated by the binary companion. 
Among the 1 1 other transiting planets with Rossiter-McLaughlin 
measurements, the only system to show a significant spin-orbit 
misalignment is XO-3, another massive and eccentric planet 
(Hebrard et al. |2008l Winn et al. |2009] l. HAT-P-2b is aUgned 
(Winn et al. l2007l LoeiUet et al. l2008l l. 

HD 80606 b is thus a new Rosetta stone in the field of plane- 
tary transits. By orbiting a bright star (V - 9), it opens opportu- 
nities for numerous follow-up studies, including: observation of 
a full photometric transit from space or multi-site campaigns to 
measure a complete spectroscopic transit sequence. 

After submission, we learned that other observations of this 
system had been obtained on the same day. Fossey et al (2009) 



and Garcia-Melendo & McCuUough (2009) independently con- 
firm the detection of the photometric transit. In addition, MEarth 
observations (D. Charbonneau, priv. comm. and oklo.org) held 
in Arizona show a flat lightcurve of HD 80606, which limits the 
transit duration to less than 12 hours, hence reinforcing evidence 
of a spin-orbit misalignment (solutions of red and blue models 
in Fig. [1] are rejected). The grazing eclipse configuration would 
also result in a slightly larger planetary radius. Further analyses 
will follow in a forthcoming paper. 
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Fig. 4. Radius of transiting gas giant planets as a function of 
the equilibrium temperature {T^q ~ r»(/?»/fl)'^^(l - e^)"'^^). 
HD80606 is circled. Its position reinforces the correlation be- 
tween incident flux and radius. All transiting gas giants are in- 
cluded above 0.4 Mjup. The red circles show the planets with 
mass higher than 2 Mjup. They follow the same tendency as 
Jupiter-like planets (black squares). 
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